Fundamental thermal fluctuations in microspheres 



o 
o 

(N 

Ph. 
< 



o 



X3 

o 



> 

o . 

o 
m 
o 



I 

o 
o 



X 
S3 



M.L. Gorodetsky and I.S. Grudinin 
Moscow State University, Faculty of Physics, Moscow, Russi^ 

We present theoretical analysis and the results of measurements of fundamental thermorefrac- 
tive fluctuations in microspheres. Experimentally measured noise spectra are consistent with the 
theoretical model. 
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INTRODUCTION 



The development of modern technology in many fields 
leads to further miniaturization of components. This 
makes it necessary to take into account certain funda- 
mental physical limitations. An example of such limita- 
tions are thermodynamical fluctuations of temperature 
in a small volume. These fluctuations are transformed to 
wideband noise in the output channel due to the temper- 
ature dependence of device parameters. The same lim- 
itations also frequently appear in experimental physics 
in macroscopic high-precision measurements. As it was 
recently shown, fundamental fluctuations with the same 
origin limit the sensitivity of gravitational wave antenna 
(international project LIGO'l"|), where thermal expan- 
sion, thermal dependence of refractive index and Young's 
modulus give rise to different types of noise Q, H, 
Thermorefractive fluctuations lead to phase noise in long 
fibers HQ , which were used to observe the effect for the 
first time 



Microspheres ^ are relatively novel type of optical 
resonators uniquely combining small size (from tens to 
thousands of micrometers) and high quality-factor up to 
Q ~ 10^° for the so-called "whispering-gallery modes" 
(WGMs). The small size of the effective volume, occu- 
pied by the e.m. field of the mode, makes for low thresh- 
olds of nonlinear effects such as bistability and oscilla- 
tory instability [lol |. which are preconditioned by Kerr's 
and thermal effects. The small volume makes it pos- 
sible for such resonators to be used as a tool for de- 
tection and measurement of the thermorefractive noise 
Apart from its importance in the LIGO project, 
the measurement of thermorefractive noise can serve as 
an innovative experimental examination of the theory of 
microscopic fiuctuations of temperature. Thermorefrac- 
tive noise should also be taken into account in possible 
applications of whispering-gallery modes resonators such 
as diode laser stabilization. 



II. THERMOREFRACTIVE NOISE IN 
MICRORESONATORS 

One can understand the effect from the well known 
thermodynamical equation for the variance of tempera- 
ture fluctuations u in the volume V: 



(1) 



where T is the temperature of the heat-bath, k is the 
Boltzmann constant, p is density, and C is specific heat 
capacity. By substituting in parameters for fused silica: 
p = 2.2 g/cm^, C = 6.7 X 10^ erg/(g • K) and effective 
volume of field of the most localized WGM in the mi- 
crosphere of radius i? ~ 50 pm with Ve/f — 10~^ cm'^, 
we obtain the value of the standard deviation of tem- 
perature a/ (u^) ~ 30 pK. These temperature fluctua- 
tions, combined with the coefficient of thermal refraction 
dn/dT = 1.45 x 10"^ K'^ lead to the effect of relative 
eigen frequency fluctuations Slu/uj ~ 3 x 10~^°. The es- 
timate is comparable with the bandwidth of resonances 
achievable in microspheres. We do not consider thermoe- 
lastic noise in microspheres, since the coefficient of ther- 
mal expansion a = 5.5 x 10~^ is sufficiently smaller 
in fused silica than dn/dT. In subsequent sections we will 
present more rigorous analysis aimed at finding spectral 
properties of this noise, taking into account peculiar field 
distribution of WGM. 

Variations of the refractive index n in the dielectric 
cavity perturb its resonant frequencies. The perturbed 
wave-equation has the form: 



^2 

Ai? + (e° + 2n5n)—E = 0, 



(2) 



where E is the electric field strength in the cavity, e*^ — v? 
is permittivity, and 5n — Jyu is the variation of refrac- 
tive index due to fluctuations of temperature u. li Eq 
is the ortho-normalized field distribution of an eigen- 
mode of the unperturbed cavity (/ EiE*dr = Sij) and 
uj = ujQ + Slo is the frequency shift then after multiplica- 
tion of this vector equation on complex conjugated vector 
Eq and integration over the whole volume, neglecting the 
second order terms we obtain: 
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where u is the temperature deviation, averaged over the 
mode volume. 



III. POWER SPECTRAL DENSITY OF 
THERMAL FLUCTUATIONS 

To calculate the effect of fluctuations of temperature 
the method of fluctuational thermal sources F(r, t) may 
be used HQ: 

fill 

^-DAu = F{r,t), (4) 

where D = X* / (pC) is thermal diffusivity and A* is 
thermal conductivity (A* = 1.4 x 10^ erg/(cm • s • K), 
D = 9.5 X 10~'^cm^/s for fused silica). This approach is 
analogous to the Langevin approach, which uses fluctua- 
tional forces in the equations of dynamics. It was shown 
before that if proper normalization of the sources is used: 



where 



= {Fir,t)Fir^,t')) 
pC 



(5) 



this approach leads to the correct results, which satisfy 
the Fluctuation-Dissipation theorem (FDT). In particu- 
lar, it was shown that thermoelastic noise is associated 
through FDT with thermoelastic damping. It is also pos- 
sible to show that thermorefractive noise is connected 
through FDT with electrocaloric losses ^2 ■ 

Thermodynamical fluctuations of temperature, aver- 
aged over the mode volume, may be calculated as: 

u = J uif,t)\Eoif)\^dr, (6) 
where J \EQ\'^df— 1 and 

^ ' Dp^+tn (2^)4' 

is the general solution of Spectral correlations of 

fluctuational forces satisfy the following condition: 

B^^ = (F0',n')F*0,n)) 

, 2kT^D 



pC 



-f3H{(3- (3')S{n-n'). (8) 



We may now calculate the following averaged value 
— {u{t)u{t + t)) (correlation function of tempera- 
ture fluctuations averaged over the mode volume) and 
so, from Wiener-Hinchin, theorem the one-sided (hence 
additional factor 2) power spectral density Su{il.) of fluc- 
tuations of temperature: 

4kT^ D 
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G0) = / lEol^e-'^f^dr 



(10) 



is normalized spatial spectrum of the energy distribution 
in the resonator. To verify this useful expression we may 
integrate it over all frequencies: 



dn _ kT^ 
'2^ ~ 
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(11) 



Comparing the final expression with JQ) we deduce that 



ye7f-J\Gm'(^ = J\m\' dr 



(12) 



The same general expression for the effective volume of 
the mode in microresonator appears in the analysis of 
nonlinearity IQj and scattering in microspheres. 

It is important to note that the above expressions were 
obtained ignoring boundary conditions. If the field, is 
concentrated near thermally isolated surface, as in the 
case of microsphere, © should be modified, substitut- 
ing in © [S{(3± - /?!) + 5{f3± + ^1)] (5(/3|| - /3[|) instead 

of (5(/? — /?'), where (3± is the component of wave-vector 
of fluctuations normal to the surface and are compo- 
nents, parallel to it. 

In the analysis above, the medium was considered 
infinite (the volume of field localization is significantly 
smaller than the size of a device). Discrete spectrum of 
thermal waves should be considered for more accurate 
calculations, especially at very low frequencies: 
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where Gi, arc the coefficients of field's intensity decom- 
position into normal thermal waves ^[^(r) of the finite 
medium. 
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Now, as before, we may calculate the correlation func- 
tion of relative frequency fluctuations and power spectral 
density: 



Sum 



AkT'^D ^ 



pC 



(17) 



Applying the obtained expressions to the microsphere 
leads to the following expression for the power spectral 
density of relative frequency fluctuations: 

(0\ - ^ fdnY_ KT^Ve fdn 



1 



1 



VT~feVd^(i + (r!Tf,)3/4)2' 



(18) 



where to, 6, d, Tf, are parameters, determined by the 
mode in microsphere (see Appendix A). This expres- 
sion was experimentally verified in this work. 

It is essential to note that in a finite body S'a(ri), con- 
trary to (|18|l . remains limited at zero frequency. How- 
ever, instead of the sum (|A16|I . we used the former ex- 
pression H18|l , obtained from the continuous spectrum, to 
evaluate the results of the measurements. There were im- 
portant reasons for doing this. First of all, too many ap- 
proximations were used when obtaining the sum. These 
are crude approximations especially for lower frequen- 
cies. However, it is shown in the Appendix and con- 
firmed by numerical calculations that the sum ljA16p and 
integral ljA5|l lead to the same asymptotic dependence 
at high frequencies. Secondly, microspheres during fab- 
rication are formed on short fused silica stem which con- 
ducts heat and couples microspheres to heat-bath. Ther- 
mal exchange with the atmosphere due to convection is 
also not negligible. Therefore, microspheres are not ther- 
mally isolated. Thirdly, other parasitic noise effects dom- 
inate at low frequencies. Throughout the experiment, 
we only considered this phenomenon at frequencies of 
n > 1/tr = D/R^ - 400c-i for {R = bQpm). 



IV. EXPERIMENTAL SETUP 

The idea of measurement of thermorefractive noise in 
optical microspheres is quite simple. If one tunes the 
measuring laser's frequency on the slope of resonance 
curve of WGM, that is on the range, where amplitude of 
the signal being measured sharply depends on frequency, 
then the trembling of resonator's eigen frequency trans- 
form into fluctuations of intensity of output radiation. 
These intensity fluctuations are recorded and further pro- 
cessed. To measure thermorefractive noise we used small 
2R ~ 80 260/iTO fused silica microspheres with qual- 
ity factors of Q ~ 10^, manufactured with the hydrogen 
miniburner. In order to avoid degradation of Q-factor 
due to adsorption of atmospheric water P| microspheres 
were placed upon fabrication to the atmosphere of dry 
clean nitrogen into the special chamber where all mea- 
surements were conducted. 



Calibration signal 
1002Hz 



Locking modulation 
91.24KH7 



mm 



Vacuum chamber 



ADC 




FIG. 1: Scheme of experimental setup for the measurement 
of noise spectra in microspheres. 



WGM in microspheres were excited with the prism 
coupler using He-Ne laser (A = 0.63/im) with piezo- 
driven front mirror as in [l3 |. Laser frequency could be 
tuned within the range of approximately 0.8GHz. The 
output power of laser (BiOfiW) was reduced to Ab^W 
passing into the chamber. The reduction of power was 
necessary to weaken thermal and Kerr nonlinearities [inl | 
and was carried out with the neutral glass filters. Fre- 
quency of the laser was locked at the slope of resonance 
curve. This was done with the use of mini-spiral heater, 
placed nearby the resonator, which actively stabilized 
the temperature of the microsphere, hence stabilizing its 
eigen frequency. To provide the feedback loop a weak 
signal at 9l.2kHz was used to modulate laser frequency. 
This frequency was chosen to coincide with one of acous- 
tic resonances of the front mirror to facilitate laser fre- 
quency modulation. This signal is demodulated by the 
WGM and error signal was used to correct the current 
which feeds the heater. When the WGM is locked to 
the laser, the output of the microsphere comprises the 
thermorefractive amplitude noise. The output from the 
resonator was registered with Thorlabs PDA500 ampli- 
fied GaAsP photodetector, with AQdB gain in amplitude 
and fd — AbkHz bandwidth at this gain setting. 

The signal from detector was digitized using ADC 
(ICP DAS PCI-1802i computer board) with digitiza- 
tion speed of 333 kilosamples per second and 12-bits per 
sample resolution. Continuous sets of 6553600 (19.66s) 
points were recorded in each session of measurement for 
further processing in each session of measurement. To 
find the estimate of noise frequency spectrum, the record 
was subdivided into many equal time intervals. Power 
spectral densities were calculated using Fast Fourier 
Transform (FFT) algorithm for each interval and then 
averaged over the entire set of intervals within each ses- 
sion. 100 intervals with 65536 = 2^^ points in each were 
used to produce estimate of spectra with fine resolution 
of w 5Hz. To avoid the effect of frequency aliasing, the 
frequency band of the signal incoming from the detec- 
tor was limited by the Butterworth filter of the eighth 
order. At stop-band frequency of 166kHz (the Nyquist 
frequency in our case) this filter attenuates the signal to 
-48dB level. First order RC-filter with time constant of 
T ~ 2s was used to get rid of the constant offset. 
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A. Calibration of spectra 

In order to obtain absolute values of spectral densities 
for measured spectra, we have developed and used the 
method of calibration based on laser frequency modula- 
tion. In this method, weak sinusoidal voltage of known 
amplitude and frequency is admixed to piezo-actuator 
of laser's mirror, thereby producing weak laser's fre- 
quency modulation. Since the measured spectra repre- 
sent relative frequency fluctuation, this additional sinu- 
soidal modulation of laser frequency results in a narrow 
peak in frequency spectrum. If the amplitude of fre- 
quency modulation is known, it is possible to find the 
value of spectral density corresponding to the spike. 

To find the coefficient of transformation of voltage to 
laser frequency modulation and to measure the quality 
factors of the resonators, we performed calibration of 
piezo-actuator using following technique. If laser fre- 
quency is swept by the saw-type voltage with known 
amplitude in the vicinity of eigen frequency of the micro- 
sphere and additional phase modulation with frequency 
higher than the bandwidth of resonant curves is applied 
to the pump, one can see on the detector two addi- 
tional resonances apart from the central frequency. These 
sidebands produce frequency scale for the calibration of 
piezo-actuator and for the measurement of Q-factor from 
the width of resonant curves (see figEJ- The coefhcient 




FIG. 2: Oscillogram of frequency marks obtained with side- 
band modulation technique and saw-tooth voltage at the 
piezo-actuator as used for caUbration of actuator and quaUty- 
factor measurement. The frequency of modulation is 5AIHz 
and the corresponding Q-factor, obtained from resonant 
curve width is 7 x 10* 

of piezo-actuator has been measured and is equal to 
|| = (3,63±0.1)i^. 

If the amplitude of frequency modulation produced by 
piezo-actuator is known, the amplitude of the peak it 
produces in calculated power spectral densities can be 
used for their absolute calibration (see Appendix B). 

To produce calibration peak we applied 5~l5mV mod- 
ulation at 172.9Hz or 1002Hz (preferred frequency at 
later stages of experiment) to piezo-actuator of laser mir- 
ror, which resulted in relative frequency modulation of 
P = duj/uj = (4 ~ 11) X 10~^^. These frequencies were 



chosen to match that of bins of discrete Fourier trans- 
form of 65536 points (namely 1001.99 = 197 • i^s/65536 
and 172.93 = 34 • Fs/65536). The higher frequency of 
Ik Hz was chosen so that the calibration peak to be sep- 
arated from low frequency setup noises in the frequency 
spectrum. It is important to note that the described cali- 
bration method does not require the quality factor of the 
microsphere as well as the transformation coefficients and 
gains in electrical tracts to be known. 

After amplitude calibration is performed, additional 
correction is required to compensate for frequency filters 
in the detector by dividing the resulting spectrum by 
1 / \/l + {f / fdY — the frequency response of the detec- 
tor. Frequency response of anti-aliasing filter was also 
compensated. To do this, it was measured and approx- 
imated with the polynomial of 7-th order so that the 
compensation procedure could be carried out by point- 
to-point multiplication. 



B. Excess noises in the measurements 



In this research, the superfluous noises were repre- 
sented by the electrical interference (generally 50Hz and 
harmonics), acoustical and seismical noises, electromag- 
netic foil of devices and also by inherited amplitude and 
frequency noises of the laser. The frequency noise of the 
laser has turned out to be the most essential source of 
influence on the quantity being investigated. We have 
measured this noise together with all other technical 
noises of the setup. The same technique as for ther- 
morefractive noise measurements was used with only dif- 
ference of using relatively large microspheres with di- 
ameters of 481, 570, 508, 588, 894, 619/iTO and modes with 
large I — m numbers of the order of 100 and above. In 
these modes, thermorefractive noise appears to be small 
enough for laser noise to dominate. Measurements of the 
laser noise for all aforementioned microspheres coincided 
well enough. The microspheres operated as the frequency 
discriminators in this case. 

Laser noises calibration was carried out in the same 
way as it was done for thermorefractive noise calibration 
— with the use of calibrating peak. In measurements of 
the technical noise spectra calibrating spike at frequency 
of 172Hz was used. 

Special measures were undertaken to weaken electro- 
magnetic hindrances in the setup. We have applied 
proper signal shielding, moreover, we have shielded ADC 
board inside the computer, which resulted in the setup 
noise decrease of more than order of magnitude. The dy- 
namic range of digitization system exceeded 90dB (tak- 
ing into account averaging over many spectra). 
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C. Identification of WGMs and computer 
processing 

Upon completion of recording process, the identifica- 
tion of WGM indices was carried out. The knowledge 
of the modes' indices allows one to calculate the effec- 
tive volume, occupied by the electromagnetic field of the 
mode, as well as theoretical energy distribution within 
the volume and therefore the theoretical spectral den- 
sity. To perform the identification, filming with the dig- 
ital camera of the speckled image of the mode was car- 
ried out. The speckled image is called forth by scatter- 
ing of mode's e.m. field at the residual molecular inho- 
mogeneities on the surface of the microsphere pj. The 



V. THERMOREFRACTIVE NOISE SPECTRA 

Figure ^ presents typical high resolution spectrum of 
measured noise in a microsphere 138/im in diameter and 
the spectrum of technical noises in large 894/j,to sphere 
measured as described in llV ij| Dashed line shows theo- 
retical curve for the spectrum calculated from recognized 
value of£ — m = 4±l. Good agreement of experimental 
data with theory for the frequencies > 100-ffz is clearly 
observable. It can also be seen that the level of combined 
technical noises, which include frequency and power fluc- 
tuations of the laser and electrical noises in circuits is 
nearly 20dB lower than the observed effect, except for 
isolated spikes. Therefore, technical noises do not pre- 
vent the experimental test of thermorefractive model. 




FIG. 3: The resonator 570/im in diameter with Z — m ~ 13 
mode (magnification 88x). 

approximation formulas were derived, linking the micro- 
sphere's radius with i index and the angular half-width 
of the WGM's belt with £ and £ — m numbers. (See Ap- 
pendix C). 

Upon recording, images were magnified and processed 
with the made-on-purpose computer program of graph- 
ical identification of the WGMs. The belt's width {£ — m 
value), simulating sphere size, slope and tilt of the belt 
could be controlled interactively, so that the theoretical 
framework projection of the mode could be made coin- 
ciding with the photographic projection of the real mode. 
This allowed one to obtain the parameters of the mode. 
The error of estimate of the £ — m value is proportional 
to y^£ — m + 1/2 i.e. the width of the mode belt, but 
is strongly determined by the quality of the pictures ob- 
tained. Modes with £ — m < 3 can usually be precisely 
recognized while for 3 < £ — m < 6 the error in obtaining 
the indices approaches 1. Index q of the mode can be in 
principle found from the optimal coupling angle of inci- 
dence of the pump beam However, this possibility 
has not been realized in the present work. We just tried 
to excite and analyze the modes with smaller q index 
which are usually characterized by higher Q and better 
coupling. 




e' o' t' 



FIG. 4: Thermorefractive noise in microsphere (138 ± 8)fim 
in diameter, I — m — i. Lower curve represents setup and 
laser noises, straight line — theoretical estimate. 

Calibrating spikes at 1002iJz for the smaller sphere 
and at 172.93Hz for the larger one can be seen on the 
graphics. At low frequencies, peaks of technical noise at 
AC SOifz-harmonics dominate. Additional sharp spikes 
could be seen at frequencies in the order of 20kHz which 
were identified as acoustical resonances of laser mirrors, 
responding to acoustical noise produced by laboratory 
equipment. One of these resonances at 92kHz was used 
in feedback loop for active mode frequency stabilization 
(see llV A)| . These spikes, except the specially excited one 
at 92kHz, are relatively small and narrow (~ 20Hz) and 
could be seen only due to high resolution of our measure- 
ments (~ 5Hz in the range of 10-^ lO^Hz). They do not 
affect the result especially if additional averaging at high 
frequencies is applied. Figure |31 depicts results of mea- 
surement of calibrated spectra of relative frequency fluc- 
tuations in four different microspheres for six different 
modes. The sizes of microspheres and estimated mode 
parameters are given on the figure. 

Data for these graphics was obtained from the data 
presented on the previous figure using uniform in log- 
scale averaging according to the following algorithm: if 
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FIG. 5: Thermorefractive noise in 4 microspheres for 6 dif- 
ferent modes. Dashed lines — theoretical curves obtained for 
recognized modes' parameters. Larger noise corresponds to 
smaller sizes of microspheres and smaller £ — m values 



fmini^fmax I fmin)^ ^ ^ fi ^ frr 



> 



' fj — rij S /' 



I, Sj 



{,f max I f mill) 

= TT-E-S"., where 



then if 

N is the number of equidistant bins {N = 200 for 
the presented figure) in log scale, rij is the number of 
points found in bin j — 0..N — 1, fmin = WHz and 
fmax = lO^Hz. This approach allows one to estimate 
better the real spectra on higher frequencies while pre- 
serving frequency resolution at low frequencies. To make 
the figure less bulky, calibration peaks were digitally fil- 
tered out while processing. 



VI. DISCUSSION 

We conclude that the results of measurements of fre- 
quency noises in microspheres confirm the theoretical 
model of thermorefractive fluctuations. This concerns 
both frequency dependence and dependence on the mode 
parameters £ and £—m. The discrepancy is less than 30% 
in the frequency range of 3 x 10^ 10^ Hz (fi > I/th), 
where the model and approximations used are valid and 
which is free of installation noises. For frequencies above 
IkHz, differences are within the limits of calibration 
and WGM identification errors. Such a close agree- 
ment was obtained in many resonators of average size 
of 120 250/xm in diameter. In very small resonators 
at frequencies lower than IkHz, the noise is up to two 
times higher, and the dependence on frequency is more 
pronounced. In large resonators at lower frequencies the 
noise seems to be smaller. However, as it was noted, the 
theoretical model described below may be not valid when 
the frequency of noise is comparable with the inverted 
time of temperature relaxation of the whole microsphere 
via the microsphere stem {Rg ^ 20^ m). The influence 
of the stem for smaller spheres due to the ratio Rg/R is 
larger. This looks compatible to the observed systemati- 



cal deviation at low frequencies of experimental curves 
from the theoretical ones obtained for thermally non- 
isolated spheres. Other mechanisms of thermal noise, 
such as thermoelastic 3 noise and thermal fluctuations 
of effective eccentricity of microspheres leading to addi- 
tional frequency fluctuations proportional to £ — m can 
also contribute to the effect. Additionally, there are other 
mechanisms of influence, which can modify spectra, such 
as convective thermal relaxation and influence of WGM 
locking scheme. 

It is worth noting, that from the displayed spectra it is 
possible to calculate the spectra of microscopic fluctua- 
tions of temperature in the volume, occupied by the e.m. 
field of WGM. Using the known relation 



6lu 
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1 dn 
n dT 



ST, 



(19) 



one can rescale the spectra with the factor ofl/ (^f^) — 
lO^K for fused silica. For instance, the mode of the 
138/im resonator has an effective volume in the order 
of 4, 3 X 10~^cm'^, which is 4 orders less than in former 
measurements of thermally induced phase noise in long 
fibers 0. 



APPENDIX A: CALCULATION OF 
THERMOREFRACTIVE NOISE IN 
MICROSPHERES 

To simplify the analysis we perform calculations below 
only for fundamental whispering-gallery mode TEm, 
which has the smallest volume of localization. The re- 
sult in principle may be extended to account for TEimq 
modes with q, £ ~ m ^ 1. The field distribution of this 
mode may be approximated as follows: 



-Icos^ e/2+u 



knR 



£+ 1/2 + 1.8558(^-1-1/2)1/3 



V;i2— Ti?3/27r3/4 

]i(knr) /]^{knR) for r < R 
g-7(r-fl) for r > i? ' 

" (Al) 



1 



However, even this approximation is too complex for an- 
alytical evaluation, so below we shall use the following 
Gaussian approximation of radial dependence: 



1 



' — 5P — 



fcni?o ~^+ 1/2 + 0.71(£+ 1/2)1/3 
d ~ Rq£-^'^ 

OMRq£-^/^. 




iiiknRo) 



(A2) 



This approximation describes rather adequately the dis- 
tribution of optical energy inside the resonator and al- 
lows one to calculate Fourier integrals. Moreover, as it 
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will be found below, because of the small depth of the 
field (parameter b), radial distribution practically does 
not influence the frequency fluctuations at frequencies of 
interest. 

G(/3) - T^^rr^ / / / e S 



xe 



2T:^bdRo 

000 

z/3r(cos 9 cos i9+sin ^ sin i9 cos((/)— (^)) ,,4^ ^2 



()3ei COS a)'' (0b sin 9)^ 



e -4 e 



1 



■n(3Ra sin i9 



sin 9d9 d(j) r dr 
(A3) 



To obtain this result, we used several approximations 
while calculating the integrals. First integration was 
done over angle 0, taking into account that the mode 
distribution is narrow and hence sin(0) ~ 1. Second 
integral over <j) leads to Bessel function of the first or- 
der, which was approximated as spherical wave. Finally, 
integral over r was found considering only fast varying 
functions. Other physical conditions were also used to 
neglect small terms: {R ~ Rq) <C R, j{knR) <^ j{knRo) 
and cos(/3i?) ~ 1 — due to boundary conditions for ther- 
mal waves. In case of /3 — > 0, the last term in ljA3p should 
be equal to 1 and the final expression is incorrect due to 
used approximations for Jq. However, this case is not 
very interesting for us. To correct the situation and ob- 
tain better approximation one may formally add 1 to the 
denominator. 

Using @ and keeping in mind additional factor 2 due 
to boundary conditions we obtain: 



SkT^D f f3'^\G{f3)\^ d(3 2kT^D 



pC J + (27r)3 pCir^R 



(f)dcoE,a)3 (0b Sin e)^ B^d'ddB 
e 2 e 2 „^ , (A4) 







J X 
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V2 



I3^£f^dp^ 

D2/34 + - AD^^Vn (1 + (f^Tb)3/4)2 ' 



(A5) 



(A6) 



This integral can be expressed through Lommel special 
functions, however, an approximation is used here, which 
works well for flTb < 1, where = {n/4y/^b'^/D. In- 
deed, neglecting first item in denominator we obtain 
high-frequency approximation with dependence fl^^ 
and neglecting exponential approximation for lower fre- 
quencies r2~^/2 - dependence is obtained. 



Finally: 



dnV KT'^Ve 

It) n^^n'^R^y/TF^CU 
1 1 

^1-byd^ (l + (f)T6)3/4)2' 



(A7) 



Numerical analysis which is omitted here, shows, that 
for the modes with I ^ m, the power spectral density of 
the fluctuations is proportional to — m) -t- 1. This 

dependence is also confirmed by the fact, that azimuthal 
"width" of the modes, and hence their effective volume, 
have the same dependence. 

To account for finite size of microspheres, calculations 
using p7|l should be performed. As at room temperature 
the power radiated from the surface (Stephan-Boltzmann 
law) is much lower than the heat exchange due to ther- 
mal conductivity (Fourier law) the simplified boundary 
condition is used below: 



9u(r, 9, 0, t) 



dr 



(A8) 



r=R 



^LM,N = CLM,mL{f3LNr)V'^ {cose) , (A9) 



c 



L.M,N 



2L + 1 {L-M)\ 
TTil + SoAi) {L + M)\ 



RHeLN-L{L- 



mmLNY 



G 



L.N 



Cl,o, 



N 



nbdRo 



(AlO) 



(All) 



~hiPLNr)PLicosl 



)r^dr sin 



0d0. 



Although only functions with M = lead to nonzero 
integrals, it is appropriate to note, that functions with 
M = 2£ could also be taken into account — they lead to 
the coupling between counter-propagating modes in the 
sphere and mode splitting. 

To estimate this integral for small values of cos 6* ~ 
-0 — 9 — tt/2 near equator of microsphere, where the e.m. 
field is concentrated, the following approximation for the 
Legendre polynomial may be used for 
L > (Po(cos6') = 1): 



Pi (cos ( 



X cos I (i+ ^)V'+ 4r 



(A12) 
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and for the spherical Bessel functions: 



jLiz) 
^LN — PlnR 
jhiiLN) 



— sin(z — Ltt/2), 
z 

tt{2N + L- 1) 



2 

(-1)^-1 



L = 2K. 



LN 



(A13) 



It is the crudeness of approximation for the roots of 
derivative of spherical Bessel functions that limits the 
applicability of the final sum we obtain below. In this 
way, the calculations below may be considered only as 
illustration. 



\Gl,n\ 

niL,N) 
7^(o,7V) 



2n{L,N) ^i^±Uf^^. 

(1-1/4L)2(1 + 1/2L) 

l-L(L+l)/ei^ 
7r/4 



if L>Q 

(A14) 



In order to calculate these coefficients the approximation 
Pln{R — Ro) ^ 1 was used. And finally: 



Sum 



TT^pDCR 



(A15) 



oo oo 



EE' 

K=aN=i 



^4 



-n{L,N). 



To calculate the sums in high-frequency approxima- 
tion, we can consider Tl{L,N) = 1. By making the fol- 
lowing substitutions: X = {K + N)/V2, y = {K-N)/V2 
and Til — R^ / D, we use integrals instead of sums: 



APPENDIX B: RELATION BETWEEN 
REGULAR CALIBRATION SIGNAL AND NOISE 
SPECTRUM 

To find the link between the amplitude of harmonic 
frequency change, corresponding peak in calculated spec- 
tral density and noise spectral density, we consider the 
Fourier transform for regular and chaotic processes. Let 
x{t) denotes chaotic process with correlation function of 

oo 

B{t) =< x{t)x(t + t) >= [ S^{n)e''^^^, (Bl) 

J 2n 



where S'^(ri) stands for "double-sided" spectral density 
which is symmetrical and two times less for positive fre- 
quencies than "one-sided" S{n). Let T denotes the 
length of the "window" W{t) in time domain, that is 
the duration of a set {x{t)} (or {xi} in discrete case) for 
which the Fourier transform is carried out. Now we shall 
consider the transform with the square window w{t) = 1 
for t e [-T/2, T/2]. Fourier transform for x{t) is: 



T/2 



X, = 



27rj 



J - I W{t)x{t)e~'^^^*dt, where ilj = 

-T/2 



(B2) 

In the procedure of spectrum estimation the spectra 
Xj{U) are calculated for small intervals of data and then 
averaged, which in our notations corresponds to the cal- 
culation of the quantity 



T/2 



{x{t)x{t'))e''^'^^-^'Utdt' 



-T/2 



S{lo) 



4sin^[(wj - oj)'^] duj 
(luj — uj)^ 2ti 



(B3) 
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47r4x4 



-X 



TT^pCDR V Stt 



t4 + 1724 



e-^''>Vi2R^)eri ( ^ ) dt. 



(A16) 



This integral coincides with ljA5|l obtained for infinite 
media, if erf-function, which is practically equal to unity 
for t > 1/VI, is ignored. However, spectral densities 
described by integral ljA5|l and sum (|A16p are different 
at frequencies for < 1/r/j where (|A16|I is finite for 
— !■ (the term with two sums is close to unity) . 



Calibrating peak in our case corresponds to frequency 
LUj >> T~^. Then, the core in the latter integral has 
sharp maximum at w = Uj and (|Xjp) = TS{ujj). 

Let now y{t) — Yq cos(r2i) represents a regular har- 
monic signal. The Fourier transform for such signal for 
square window will look as follows 



T/2 



Yq cos(f2i)e 



dt. 



(B4) 



-T/2 



leading to Ysj — YoT/2. Now we compare the results of 
spectrum estimation procedure (jB3|l and the results of 
Fourier transform for harmonic signal: 



r2 



Y'^T'^ 
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Here one can easily see the link between the amplitude 
of harmonic signal and corresponding peak in spectral 
density. If the harmonic signal has the dimension of 
frequency and the desired spectral density represents 
relative frequency changes, then the final relation, cor- 
responding to one-sided spectral density of calibration 
spike with square window in Fourier transform, will be 



S. 



1 



(B5) 



where /o — is optical frequency. The same calculations 
for the Hann window give: 



where 



tiq = + 



(C3) 



These formulas allow to get an estimate of ^ and £ — m 
values from the radius of the microsphere and width of 
the belt of the mode visualized due to residual surface 
scattering. 



'Sf/f 



Yo 
/o 



1 



(B6) 



APPENDIX C: FORMULAS FOR WGMS 
IDENTIFICATION 



To identify the parameters of whispering-gallery 
modes from photographs the following approximation for 
half-width of WGM belt was used: 



AO = V2(^+l/2-m)/£ 



(CI) 



To calculate £ index from the radius of microsphere we 
found by numerical approximation the following approx- 
imated formulas valid for i? » A: 



(7=1 

q = 2 

(7 = 3 

(7 = 4 



I = -0.5 + - 2.287t 



1/3 
Iq 



0.1718t;g^/-\(C2) 



I = -0,5 + tiq- 4:M7t 
I = -0.5 + tiq - 6.895t 



1/3 



Iq -1 

1/3 
Iq 

1/3 
Iq 



0.6944t 



l.SlSt 
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